A New Damping Mechanism in Non-linear Bubble Dynamics 
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Non-linear equations of radial motion of a gas bubble in a compressible viscous liquid have been 
modified considering effects of viscosity and compressibility more complete than all previous works. 
A new set of equations has been derived including new terms resulted from consideration of the 
viscosity and compressibility not only at the bubble interface, but also in the bulk of liquid. The 
new equations are two non-linear coupled equations, which can not be merged into one equation 
unlike all previously derived equations. Numerical calculations have been performed considering 
effects of heat and mass transfer at the bubble interface. The results indicate that the new terms 
exhibit an important damping role at the collapse, so that their consideration dramatically weakens 
the bubble rebounds after the collapse. Dependence of this new damping mechanism to amplitude 
and frequency of the deriving pressure has been investigated. 

PACS numbers: 47.55.Bx, 43.25.Yw, 43.25.+y, 78.60.Mq 



I. INTRODUCTION 

When a small isolated gas bubble, immersed in a liq- 
uid, experiences a high amplitude spherical sound field, 
it grows and contracts non-linearly. Description of the 
dynamics of such non-linear motion is an old challeng- 
ing problem. The complexities of the problem arise from 
the effects of heat conduction, mass diffusion, compress- 
ibility, viscosity, and energy losses involved in damped 
oscillations of the bubble. So far, many people have in- 
vestigated this problem, each concentrating on different 
aspects of its complexities. However, a rather complete 
description has not been presented yet. 

Lord Rayleigh Q] was the first who analytically solved 
the collapse of an empty bubble in an incompressible liq- 
uid. Plesset @ subsequently studied the dynamics of 
a vapor-filled bubble in a time-dependent pressure field. 
Noltingk and Neppiras Q were the first people derived 
a differential equation for the gas bubble motion under 
the influence of such pressure field. The extension of 
this equation to the bubble motion in a compressible 
liquid has been studied by many authors; Herring 0, 
Trilling H, Gilmore @], Keller and Kolodner 0, Hick- 
ling and Plesset Jahsman Flyrm Lastman 
and Wentzell [TJ, Keller and Miksis [l^]- On the other- 
hand, heat conduction effects were presented in the works 
of Hickling Fujikawa and Akumatsu [l4|, an d Yasui 
[l5l | . In addition, the works of Lofstedt et al. , Nig- 
matulin et al. [13, and Brujan are of the subsequent 
papers addressing this problem. 

Between all previous works, the most complete descrip- 
tion of the bubble dynamics equations was presented 
by Prosperetti and Lezzi [l^. They used a singular- 
perturbation method of the bubble-wall Mach number 



and derived the following one-parameter family of equa- 
tions describing the bubble motion in the first order ap- 
proximation of compressibility [l9j : 



R d (P p ) 




(1) 



where, R, C, P^, and p are the bubble radius, liq- 
uid sound speed, liquid pressure at infinity, and liquid 
density, respectively. Also, 77 is an arbitrary parameter. 
Equation Q must be supplemented by a boundary con- 
dition equation at the bubble interface to relate the liquid 
pressure, Pi, to the gas pressure inside the bubble. Like 
all previous authors, Prosperetti and Lezzi used the 
following incompressible equation for this purpose: 



Pi =Pn 



A R 2a 



(2) 
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where, P g , fi, and a are the gas pressure at the bubble in- 
terface, liquid viscosity, and surface tension, respectively. 
Most of the previously derived equations belong to this 
single parameter family of equations, corresponding to 
different values of 77. Moreover, rj — yields results in 
closest agreement with the numerical simulation of full 
partial differential equations 0] . 

Two specific approximations have been used in the 
derivation of Eq'ns. Q and The first approxima- 
tion is the negligence of the viscosity effects for the liquid 
motion around the bubble, which has been used in the 
derivation of Eq'n. 0J. In fact, Eq'n. has been 
derived from the Euler equation, in which the viscosity 
is eliminated. Note that the viscous term of Eq'n. 
has been resulted from the liquid viscosity at the bubble 
interface, but not from the bulk of liquid. 
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The second approximation is the incompressibility as- 
sumption of the liquid and the gas at the bubble inter- 
face, which has been used in the derivation of Eq'n. J3J). 
All of the effects of the liquid compressibility in the work 
of Prosperetti and Lezzi, as well as in all other previous 
works, have been resulted from the liquid motion around 
the bubble, but not from the bubble boundary condi- 
tion equation. In fact, all previous authors, on one hand 
took into account the compressibility of the liquid mo- 
tion around the bubble, but on the other hand neglected 
its consideration at the bubble interface. 

Although, the two mentioned approximations have 
been used in the derivations of the existing bubble dy- 
namics equations Qp, USEHHIEHIIflmElil 

Q [H HI [H [llE37 but the applicability of these ap- 
proximations for all times of the bubble motion needs to 
be clarified. Especially, at the end of the collapse, when 
the bubble motion is significantly compressible, these ap- 
proximations may be inapplicable. 

In this paper, by eliminating the above mentioned ap- 
proximations, we have originally modified the bubble dy- 
namics equations. A new set of equations have been de- 
rived including all effects of the viscosity and compress- 
ibility of both the liquid and the gas. These equations 
contain new terms resulted from the effects of two coef- 
ficients of viscosity of both the liquid and the gas. The 
influence of the added new terms arisen from the liquid 
viscosity has been numerically investigated. The results 
clearly indicate that the addition of the new terms con- 
siderably affects the bubble motion at the collapse time 
and during the bubble rebounds. 



II. DERIVATION OF THE BUBBLE DYNAMICS 
EQUATIONS 

To derive the bubble dynamics equations, we assume 
that the motion of the bubble interface and the infinitely 
extended surrounding liquid are always spherically sym- 
metric. Under this circumstance, the continuity and the 
momentum equations are as follows 20] : 
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where, p, u, p, and A are density, velocity, pressure, and 
divergence of the velocity, respectively. Also, u and A 
are first and second coefficients of viscosity. The two 
brackets of the right hand side of Eq'n. J3J have the 
same differential forms. Hence: 

f du du\ dp n . dA 



Dividing Eq'n. by p and integrating it with respect to 
r from the bubble interface to infinity at a fixed instant of 
time and assuming that the two coefficients of viscosity 
are constant, we get: 



dt r 2 



POO -l POO 

/ -dp + (A + 2u) I - dA, 
Jr P Jr P 



(6) 



where, the liquid velocity at the far field is assumed to 
be sufficiently small. When the liquid density does not 
strongly change, Eq'n. © can be written as: 



du R 2 P^-Pi \ + 2p , A 
r dt 2 p p 



(7) 



where, A/ and A^ are the divergence of the liquid veloc- 
ity at the bubble interface and at infinity, respectively. 
The liquid pressure at the far field, Poo, is: 



Poo(<) =P +Pa(t), 



(8) 



where, Pq and P a (t) are ambient and driving pressures, 
respectively. From Eq'n. (J3J) the divergence of the veloc- 
ity can be written as: 



A = -l^E = L^E 

p dt pc 2 dt 



(9) 



where the sound speed, c, which is assumed to be con- 
stant in the liquid, is defined as: c 2 = dp/ dp. Thus, Eq'n. 
Q) becomes: 

du, R 2 Poo - P, A + 2/i d 
„ -d-t^^^-^p—^^dt^-^- 

(10) 

The viscous term in Eq'n. (|10fl has been resulted from 
the simultaneous effects of the liquid viscosity and com- 
pressibility. This term has been neglected in the derived 
equations of the previous authors. 

To convert Eq'n. IjlOl) to a differential equation, the 
integral term must be approximated. Because of the irro- 
tationality assumption of the bubble motion the velocity 
potential, ip(r, t), can be introduced into this equation as: 

2 p p A C A dt 

(11) 

where subscript t denotes the temporal derivation. In the 
simplest approximation of ip t (R, t) , the liquid motion is 
assumed to be incompressible. So, the liquid velocity at 
any distance r will be: 



u{r) 



R R 



(12) 
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According to this approximation, (pt(R, t) can be written 
as: 



<p t (R,t) 



'1 




dt 





= RR + 2R* 



(13) 



J r=R 



Inserting Eq'n. (|13|l into Eq'n. Ullfl and eliminating the 
liquid viscous term due to incompressibility approxima- 
tion yields the well-known Rayleigh- Pies set equation: 



RR + Itf = __L__* 

2 p 



(14) 



To introduced the compressibility effects into approx- 
imation of ip t (R,t), we assume that the liquid potential 
fulfills the following acoustic equation for the spherical 
waves: 



0_ 

dt 



(rep) 



or 



0. 



(15) 



As we have shown in the following, this approximation is 
equivalent to the first order compressibility consideration 
of Ref. ^| . Double radial differentiation of Eq'n. Ill 511 
results: 



|(rA) + c| : (rA) = 0. 



(16) 



Through the definition of total derivative of the velocity 
along with Eq'n. © it can be obtained that: 



fn ^ r du j r( du a 2u2 \ , „ 



By a partial integration: 

(p t (R, t) = RR + 2R 2 - RRAi + 



°° dA , 
r——dr. 



, R dt 

Applying Eq'n. lfTo|l in Eq'n. l(T%|) yields: 

<p t (R, t) = RR + 2R 2 - RRAi + RCA t . 
Substituting Eq'n. |(T%1) into Eq'n. (fTTf) results: 



7) 



(18) 



(19) 



RR 



3.. 



:R Z = 




p 2 C 2 dt 

This equation is the modified form of Flynn equation p"fj| 
along with viscosity consideration. 

The compressibility effects can be introduced in a dif- 
ferent manner. Differentiating Eq'n. (|15fl with respect to 
t and substituting tp t (R,t) in it from Eq'n. I|ll|) yields: 



RR fdu\ 
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C \dtJ R P C 



R 



R dPg, 
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Poo -Pi _ X + 2p d 
p p 2 C 2 dt 



(Pi -Poo), (21) 



in which, the spatial and the temporal partial derivatives 
of the viscous term have been neglected due to smallness 
as well as the spatial derivative of the liquid pressure 
at infinity. Inserting the relation dp — dp/c 2 into Eq'n. 
(0 and using this equation and the momentum equation 
along with the total derivatives of the pressure and the 
velocity at the bubble interface, the four partial deriva- 
tives in Eq'n. (|21|) can be obtained as: 



du\ 
dt) R 



dp 

dt 
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dr 
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IT' (22) 



= pRR 



m 

dt ' 



1 dPi 2R 
'p~C 2 ~dT ~ ~R' 



-pR. 



(23) 



(24) 



(25) 



Note, in the derivation of the Eq'ns. I|22I25[) . the spa- 
tial and the temporal partial derivatives of the two vis- 
cous terms in the Navier-Stokes equation, Eq'n. 10}, have 
been neglected. Inserting Eq'ns. (|22I25() into Eq'n. J5TJ 
and retaining only the terms up to the order of R/C re- 
sults the modified form of Herring- Trilling equation 0,0 
along with viscosity effects: 



Pi - Poo A + 2/x d 
p P 2 C 2 dt 



(Pi - Poo) ■ 



(26) 



It is straight forward to obtain the modified form of 
Eq'n. iJTJ by the addition of Eq'n. 126|) and the product 
of (1 - rj)R/C with Eq'n. as: 



l-(, + l)§J^ + (l-(3, + l)Aj \k 2 = 

pcdt {Pi - p - ) + [ 1+ ^-^c)—r~ 



A + 2p d 
p 2 C 2 Jt 



(Pi -Poo). (27) 



The new Eq'n. l(77l) provides a suitable description of 
the simultaneous effects of liquid compressibility and liq- 
uid viscosity in the bubble motion. The added viscous 
new term of this equation arises from the liquid com- 
pressibility. 
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To complete the argument, it is necessary to modify 
the boundary condition Eq'n. (J2J with the compressibil- 
ity effects. The radial component of the stress tensor 
is: 



T rr = — p + AV • u + 2p I — 

or 



(28) 



Inserting the velocity divergence from Eq'n. into this 
equation yields: 



, , „ x / 9u 2u\ uu 
\or r J r 

= - P +(X + 2p)A-^. 



(29) 



The boundary continuity requirement at the bubble in- 
terface is: 



T rr (liquid) \r= T rr {gas) \r +2 



R 



(30) 



Substituting Eq'ns. |2HI> and © into Eq'n. (jSOJ leads 



to: 



uR 

Pi +4— 
R 



X + 2p dPi 
pC 2 dt 

Xg + 2/Ig dp g 



P„ 



R 



Pa 



dt 



4 



(31) 



where, p g , p g , and A g are the gas density, the first and 
the second coefficients of viscosity of the gas, respectively. 
The new Eq'n (|31|l provides the most complete boundary 
condition equation at the bubble interface, which con- 
tains all effects of compressibility and viscosity of both 
the liquid and the gas. Here, we concentrate on the effects 
of liquid viscous new term of Eq'n H31[) and eliminate the 



Eq'n. l]31[l becomes: 



Pi + 



A + 2/i dPi 
pC 2 dt 



pR 
~R 



-2°-. 
R 



(32) 



Comparison of Eq'ns. J2J) and (|32|l indicates the exis- 
tence of a new viscous term in Eq'n. 132|) due to the 
liquid compressibility. The set of Eq'ns. (|2*7jl and J3^|) 
present the most complete form of the bubble dynamics 
equations containing effects of the liquid compressibility 
and viscosity not only at the bubble interface, but also 
in the bulk of liquid. In fact, these equations account for 
the viscosity of a compressible liquid. While, all previous 
equations accounted for the viscosity of an incompress- 
ible liquid and compressibility, separately. 

To generalize the argument, Eq'ns. I|27|l and l|32|) 
are expressed in dimensionless forms. The dimensionless 
variables of this problem are defined as: 



R* = 



R 

Rn 



R 



tC 



R* = q, r = — , p; = 
p 

p* = Jj_ P 
9 P c 21 00 P c 21 



p, 



i 

Ro' ' l p~C^' 
p 



(33) 



where, Rq is ambient radius of the bubble. Substituting 
the dimensionless variables into Eq'ns. I|27l) and the 
dimensionless equations are obtained as: 

(l - (r, + l)i?*) R*R* + | (l - |(3r7 + 1)B*\ R* 2 



l + {l-ri)R*){Pt-Pl-P^) 



P* = p* - - 2— - (A* + 2^*)^. 

I g R * i?* dt* 



(34) 



(35) 



The quantities a*, p*, and A* are dimensionless surface 
tension and dimensionless liquid viscosity coefficients, 
which are defined as: a* = a/pRoC 2 , p* = p/pRoC, 
and A* = X/pRgC. These dimensionless numbers, which 
are basically inverse of Weber Number and inverse of 
Reynolds Number, characterize significance of the surface 
tension and the liquid viscosity in the bubble dynamics. 



III. BUBBLE INTERIOR EVOLUTION 

To quantify effects of the new viscous terms on the bub- 
ble dynamics, evolution of the gas pressure at the bubble 
interface, P g , must be specified. It can be determined 
from simultaneous solution of the conservation equations 
for the bubble interior and the bubble radius equations 
|H 0|H IHIH |2| . Also > heat conduction and mass 
exchange between the bubble and the surrounding liquid 
affect the bubble evolution. These complexities were con- 
sidered in a complete gas dynamics simulation by Storey 
and Szeri 29]. 

On the other hand, strong spatial inhomogeneities in- 
side the bubble are not remarkably revealed, unless at the 
end of an intense collapse [26l |27|. Therefore, the uni- 
formity assumption for the bubble interior seems to be 
useful and provides many features of the bubble motion 
[2ll 1221 ] . Using this assumption Toegel et al. presented an 
ODE model , m which effects of heat and mass trans- 
fer at the bubble interface have been considered. This 
model accurately describes various experimental phase 
diagrams |3l| and provides a good agreement with the 
complete direct numerical simulation of Storey and Szeri 

EH- 

In this paper, we used the Toegel et a/.'s model |3(| 
for specifying the bubble interior evolution. We describe 
an argon bubble in water under the conditions of Sin- 
gle Bubble Sonoluminescence |2ll l22j . The gas evolution 
model can be summarized as follows. 

The gas pressure is modeled by the van der Waals equa- 
tion of state: 



Pn 



N tot kT 

V - NtatB 



(36) 
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where, N to t = Nat + Nh 2 o is total number of parti- 
cles inside the bubble. The covolume constant B = 
5.1 x 10~ 29 m 3 is assumed to be equal for both water 
vapor and argon |3l| . The value of N to t changes with 
time because of evaporation and condensation of the wa- 
ter vapor molecules at the bubble interface. The rate of 
change can be modeled as [30| : 



tot ■ 



-AttR 2 D 



n 0f 



Id 




l d = mini , /— r-,— 



(37) 



where, n H ^ and no H2Q are the instantaneous and equi- 
librium concentration of the water vapor molecules, re- 
spectively. The diffusion coefficient, D, is given by |32| : 
D = D Q (n /n tot ), where D = 23.55 x 10~ 6 m 2 / s and 
iiq = 2.5 x 10 25 m 3 . The quantit y L is thickness 
of diffusive boundary layer (see Ref. [3(| for more de- 
tails). The equilibrium concentration of the water va- 
por molecules is given by the number density of the 
saturated vapor pressure at the ambient liquid temper- 
ature To; no H20 — P v {To)/kTo. In our calculations: 
n 0ll2O = 5.9 x 10 23 m" 3 . 

Similar to Eq'n. (|37|) . the heat exchange at the bubble 
interface can be approximate by |30j ]: 



Q = AttR 2 k 



Tq-T 
I 



Hh 



ill 




(38) 



where, T and To are the gas temperature and the am- 
bient liquid temperature, respectively. Also, n — 17.9 x 
1CP 3 W/mK |32( is thermal conductivity coefficient of 
the gas content, l t h is thickness of thermal boundary layer 
and x is thermal diffusivity coefficient, which is given by: 
X = Tc/cp, with c p — (^TiAr + ^nH 2 o)k the constant pres- 
sure heat capacity per unit volume of the gas. 

Applying the basic energy equation for the bubble con- 
tent evolution results the rate of change of the gas tem- 
perature |30L l3l| : 



C v = -N Ar k + 



4T -3T-T 



Bj/T 



kNtot 

c v 

(39) 
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(Oi/T) 



2 p 0i/T 



( e e,/T _ x y 



kN 



H 2 0, 



(40) 



where, three different values for 9i correspond to char- 
acteristic vibration temperatures of H2O: d\ = 2295K, 
2 = 5255K, and 63 = 5400K jp. In this paper for sim- 
plicity we neglected the effects of chemical reactions |31j , 
which arc important only at an extremely short time in- 
terval when the bubble temperature is high enough (more 

than 5000K). 

Equations (|36I4U|) along with the bubble dynamics 
equations are the set of equations, which totally de- 
scribe the evolution of the bubble characteristics. Under 





FIG. 1: Time variations of the bubble radius, according to the 
new (solid) and the old (dashed) sets of equations. Graphs (a) 
and (b) shows the bubble evolution in a complete period and 
during the bubble rebounds, respectively. The equilibrium 
radius is Ro = 4.0 fim and the deriving pressure amplitude is 
P a = 1.4 atm 



these circumstances, time variations of the bubble prop- 
erties have been numerically calculated for both the new 
(Eq'ns. I34I35|I and the old (Eq'ns. I1I2H bubble dynamics 
equations (for 77 = 0). 



IV. NUMERICAL ANALYSIS 

The calculations were carried out for a periodic driving 
pressure; P a (t) — P a sm(ojt), with u = 2ir x 26.5 kHz. 
The constants and parameters of the bubble dynam- 
ics equations were set for the water at room temper- 
ature, To = 293.15 K, and atmospheric ambient pres- 
sure, P = 1.0 atm; p = 998.0 kg/m 3 , C = 1483.0 m/s, 
fj> = 1.01 x 10~ 3 kg/ms, a = 0.0707 kg S - 2 34]. The sec- 
ond coefficient of viscosity of water at room temperature 
was set to be A = 3.43 x 10" 3 kg/ms pj . 

Figures (1-3) illustrate the results of our calculations 
for P a = 1.4 atm and Rq — 4.0 yon. Similar values for 
these parameters have been reported in recent experi- 
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FIG. 2: Time variations of the number of H2O particles inside 
the bubble according to the new (solid) and the old (dashed) 
sets of equations. Graphs (a) and (b) show the evolution in 
a complete period and during the bubble rebounds, respec- 
tively. The parameters and constants are the same as Fig. 
(I)- 



mental works of Ketterling and Apfel j3g , Simon et al. 
[33, and Vazquez et al. j^j. Figure (1) shows the varia- 
tions of the bubble radius for the new and the old bubble 
dynamics equations. It is observed that the addition of 
the new viscous terms considerably changes the bubble 
evolution after the collapse. The bubble motion is re- 
markably compressible during the collapse. Therefore, 
the new viscous terms, which have been arisen from the 
liquid compressibility, are important in this time interval. 
These terms exhibit a damping role and their considera- 
tion reduces the amplitude of the bubble rebounds. Also, 
the period of the rebounds decreases with the addition of 
the new terms. Details of our calculations at the end of 
the collapse show that the minimum radius for the new 
case is about 4% greater than that of the old one. 

The difference between the two cases for the evolution 
of the number of H2O particles has been illustrated in 
Fig. (2). As this figure indicates, a large amount of wa- 
ter vapor evaporates into the bubble during the expan- 
sion. Indeed, at maximum radius, most of the bubble 
content is water vapor. During the collapse, the water 



FIG. 3: The gas temperature evolution during the bubble 
rebounds (a) and at the end of the collapse (b) for the the new 
(solid) and the old (dashed) cases for the same parameters and 
constants as Fig. (1). 



vapor molecules rapidly condense to the liquid. It is seen 
that the remarkable difference between the two cases also 
appears on the variations of the number of H2O particles 
after the collapse. Note that, the difference gradually 
disappears as the bubble rebounds weaken. 

In Fig. (3), the gas temperature evolution during the 
bubble rebounds and at the end of the collapse have 
been demonstrated. Damping feature of the new term 
is clearly observed by considerable decrease of the peak 
temperature at the collapse time as well as remarkable re- 
duce of the secondary peak temperatures. Also, the pulse 
width of the main peak temperature decreases with the 
addition of the new terms. 

Figures (4) and (5) illustrate dependence of the bubble 
properties at the end of the collapse to the variation of 
the amplitude and the frequency of driving pressure for 
the two cases. The different values of P a , corresponding 
to a specific values of Rq, can be experimentally obtained 
by var ying the concentration of the dissolved gas in the 
liquid |2ll l36l |37], H3- Also, the frequency range in the 
calculations of Fig. (5) is in the range of the experimental 
results [2l|. 

Figure (4) shows the variation of the peak tempera- 
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FIG. 4: The bubble characteristics at the time of collapse FIG. 5: The dependence of the bubble characteristics at the 

as a function of the driving pressure amplitude for the new end of the collapse to the frequency of driving pressure for 

(solid) and the old (dashed) bubble dynamics equations; peak the new (solid) and old (dashed) bubble dynamics equations; 

temperature (a), peak pressure (b), and minimum radius (c). peak temperature (a), peak pressure (b), and minimum radius 

The equilibrium radius was fixed (Ro = 4.5 fim). Other con- (c). The equilibrium radius and the pressure amplitude were 

stants are the same as Figs. (1-3). fixed (Ro = 4.5 /J,m and P a = 1.35 atm). Other constants are 

the same as Figs. (1-3). 



ture, the peak pressure, and the minimum radius as a 
function of the driving pressure amplitude. The ambient 
radius was fixed (Ro = 4.5 /im). The damping feature of 
the new viscous terms is clearly observed in this figure. 
For the peak temperature in Fig. 4(a), the maximum 
relative difference appears around P a = 1.25 atm (about 
35%). The difference between the two cases decreases for 



higher driving pressures. However, for the peak pressure 
in Fig. 4(b), the difference increases with the increase 
of the amplitude. The values of R m in for the two cases 
are similar in high driving pressures because of the ef- 
fects of the excluded volume, which prevent the bubble 
compression |33| . 
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Figure (5) represents the dependence of the bubble 
characteristics at the end of the collapse to the varia- 
tions of the driving pressure frequency. The ambient ra- 
dius and the amplitude were fixed in these calculations 
(Rq = 4.5 fim and P a = 1.35 atm). As Fig. (4), the 
damping feature of the new terms is also seen here. The 
difference between the two cases for the peak temperature 
and the minimum radius increases when the frequency is 
increased. While, for the peak pressure the difference 
reduces with the increase of the frequency. 

A major deficiency of the old bubble dynamics equa- 
tions is that for strongly driven bubbles, such as sonolu- 
minescence bubbles, large amplitude rebounds are pro- 
duced after the collapse, so that they often last until 
the next acoustic cycle of the periodic driving pressure. 
This is in contrast with the experimental results, which 
show rapidly damped rebounds 0, |3^| . By introducing 
a dam ping term arisen from the gas compressibility, Moss 
et. al |39( provided a typical solution for this problem. 
The effects of the suggested term by Moss et. al is very 
similar to the damping effects of the new terms in this 
paper, (compare Fig. 1(b) with Figs. (3) and (4) of Ref. 
|39|). It seems the damping feature of the the new terms 
is a better way for solving the mentioned problem. The 
reason is that Eq'ns. (|34|l and 1)35(1 have been derived di- 
rectly from the basic equations of fluid mechanics, on the 
contrary to Eq'n. (3.2) of Ref. [3!|, which was derived 
by an approximate method. 



V. CONCLUSIONS 

The equations of the bubble motion in a compressible 
viscous liquid were newly derived from the full Navier- 
Stokes equations. These equations contain two similar 
new terms resulted from the simultaneous effects of the 
liquid viscosity and compressibility. These new terms 
have a considerable damping role at the collapse, when 
the bubble motion is significantly compressible. This 
new damping mechanism dramatically changes the bub- 
ble properties at the end of the collapse and during the 
bubble rebounds. 

The results of this work indicate that, the neglect of 
the new terms in the previous works is not reasonable 
for the collapse time and the new effects should be con- 
sidered for the prediction of the quantities related to the 
collapse, e.g., the value of light emission by a single sono- 
luminescing bubble as well as the bubble stability limits. 
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